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Abstract. The electromagnetic and the muonic longitudinal profile at production enclosure important informa- 
tion about the primary particle and the hadronic interactions that rule the shower development. In fact, these 
two profiles provide two different insights of the shower: the electromagnetic component gives a measurement 
of the energy and the strength of the neutral pion channel; while the muonic profile, being intimately related 
with the charged mesons decays, can be used as a direct probe for the high energy hadronic interactions. 
In this work we explore the interplay between the electromagnetic and muonic profiles, by analysing their phe- 
nomenologic behaviour for different primary masses and energies, zenith angles, and also different high energy 
hadronic interaction models. We have found that the muonic longitudinal profile at production displays univer- 
sal features similar to what is known for the electromagnetic one. Moreover, we show that both profiles have 
new primary mass composition variables which are fairly independent of the high energy hadronic interaction 
model. 

Finally we discuss how the information in the electromagnetic and the muonic longitudinal profile can be 
combined to break the degeneracy between the primary mass composition and the high energy hadronic physics. 
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1 Introduction 

The Ultra High Energy Cosmic Rays (UHECRs) are the 
most energetic particles known. And yet its origin and 
composition remains a mystery although they were dis- 
covered half a century ago. The reason for this resides 
on the difficulty of their detection. Since their flux at the 
highest energies is very low the only way to detect them 
is by observing the huge cascades of secondary particles, 
known as Extensive Air Showers (EAS), that are created 
by the interaction of these cosmic rays with the atmo- 
sphere molecules. 

The EAS can be detected by sampling the charged par- 
ticles that arrive at the ground or, in moonless nights, the 
shower longitudinal profile can be sampled through the 
collection of the fluorescence light produced as the shower 
develops into the atmosphere. The latter measurement, al- 
though with a duty cycle of around 10%, provides a di- 
rect measurement of the electromagnetic shower compo- 
nent development. 

While the aiTival direction of the UHECRs can be 
easily obtained by the experimental techniques aforemen- 
tioned, the determination of the nature of the primary par- 
ticle is much more difficult. In fact, the shower observ- 
ables connected to the type of particle are also sensitive 
to the physical interactions that rules the shower develop- 
ment. Moreover, the hadronic interactions at high ener- 
gies are described through phenomenological models that 
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are fitted to the available accelerator data and extrapolated 
several orders of magnitude to the UHECRs energies. This 
implies that the determination of the UHECRs mass com- 
position is linked to the understanding of the shower phys- 
ical mechanisms, in particular to the hadronic interactions 
descriptions. 

Muons are a fundamental tool to assess the high en- 
ergy hadronic interaction models as they are produced es- 
sentially from charged mesons (secondary products from 
hadronic interactions). Moreover, a large fraction of the 
muons can reach the ground due to a low interaction cross- 
section and a relatively large decay time. Recent analyses, 
that use the muon arrival time, are able to reconstruct the 
muon longitudinal profile at production at ultra high ener- 
gies HI], providing an additional insight of the shower 

The phenomenological behaviour of the muonic longi- 
tudinal profile at production and its relation with the elec- 
tromagnetic shower profile is the main focus of this work. 
Moreover, we show, at the end of this paper, how several 
shower observables may be combined to disentangle the 
mass composition from the hadronic interaction models. 

The paper is organised as follow: in section |2] the 
muonic longitudinal profile at production is presented and 
the main observables, namely Xi„„^, are discussed. When- 
ever relevant the features of this profile will be compared 
to the ones found in the electromagnetic profile; In sec- 
tion [3] we show that the muonic profile shape exhibits an 
universal behaviour and its coiTelation with the shape of 
the electromagnetic profile is assessed; In section (sec. 3} 
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Figure 1. Muon production shower profiles as a function of 
depth (X), for proton (red) and iron (blue) primaries at £ = 10" 
eV. In this picture 100 showers for each primary are shown, both 
generated with QGS Jet-II.03 and with 6 = 40° . 



we show how hadronic interaction models can be distin- 
guish independently of primary mass composition scenar- 
ios, through the combination of different shower observ- 
ables. Finally, some conclusions are drawn in section |5] 

2 Longitudinal shower profiles 

To study the truJ^ Muonic Production Depth profile (on- 
wards referred as MPD profile) we have used Conex ^ 
simulations (for more details we refer the reader to [fl). 

The MPeI profile for proton induced showers (in red) 
and iron primaries (in blue) at £ = lO'^ eV is shown in 
FigU Compared to proton showers, iron initiated showers 
have more muons, which is readily seen by looking at the 
maximum of the profil^E N^ax- Moreover, the height of 
the maximum fluctuates more for proton induced showers 
than for iron showers. This is a distinct feature from the 
electromagnetic profile that, having its N^^'^x related to the 
shower energy, has essentially the same maximum value 
for both proton and iron shower, whereas the fluctuation 
level is of the order of 5%. 

The second feature worth noticing is the depth of the 
shower maximum, X^,„jt, and the corresponding distribu- 
tion which can be seen in more detail in Fig.|2] Here one 
can see that the behaviour of this observable follows the 
one found in the electromagnetic profile. The X^„„ ^ of pro- 
ton induced showers are deeper and fluctuate more than in 
iron showers. This observable is sensitive to both primary 
mass composition and hadronic interaction physics, being 
the difference between proton and iron bigger than the dif- 
ference for different hadronic models (QGS Jet-II.03 



the experiments measure the apparent-MPD, which depends on the 
point of observation, as the true distribution has to be convoluted with 
the muons propagation. By its turn, the muon propagation depends of the 
distance travelled by the muons, their energy, among others. The trans- 
formation of the apparent-MPD into the fme-MPD is described in Ql- 

^note that in Conex there is an energy threshold for muons of 1 GeV. 

^in this paper the index fi and e.m. will be used to address to vari- 
ables related with the MPD profile and the electromagnetic profile, re- 
spectively. 




Figure 2. X'^„„^. distributions of muon production longitudinal 
profiles for different primaries and hadronic interaction models 
- the red lines are for proton and the blue ones for iron pri- 
maries. The distributions generated with QGSJet-II are shown 
as full lines while the ones generated with EPOS are displayed 
with dashed lines. 



results are displayed as full lines and EPOS 1.99 as 
dashed lines). 

An interesting correlation to look at is the one between 
the X^ax and X,^'^;., in an event-by-event basis. This profile 
is shown in Fig.|3]for proton (in red) and iron (in blue) in- 
duced showers, both generated with QGSJet-II. The elec- 
tromagnetic maximum is reached around 200 g cm^ later 
than the muonic profile maximum. On average, there is 
a correlation between the electromagnetic and the muonic 
Xf„ax- It depends slightly on the primary mass composition 
and it was checked to be nearly independent of hadronic 
interaction models. 

The difference between electromagnetic and muonic 
Xmax, shown as an inset plot of Fig. [3] also changes be- 
tween primaries - by around 30 g cm^ in the average sep- 
aration - which means that although related, the two pro- 
files give some independent information. This quantity, 
^m'ax~^max, IS Very interesting as it is by construction in- 
dependent of the first interaction point, which means that 
it should not be affected by sudden changes in the primary 
cross-section as suggested in fl. 

3 Universality of the longitudinal profile 
shape 

In Fig.|4]the profiles shown in Fig.[T]are expressed in X' = 
X - X,„ax and A^' = N/Nmax- In these variables one can 
see that the obtained shape is rather universal, similarly to 
what happens to the electromagnetic profile. 

Taking advantage of this universality it can be useful to 
use the average shape in order to determine the two main 
parameters from a fit with a reduced set of data. More- 
over, there can be extra information on the average shape 
similar to what happens in the electromagnetic profile [9]. 
In fact, it was found that iron showers have a broader 
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Figure 3. Relation between muon production and electro- 
magnetic X'lj '^^, for E = 10" eV showers. The difference between 
the two X„,aj are shown in the inset plots. The dependence on pri- 
mary mass is shown (proton is the red (full) line and iron the blue 
(dashed) line, both generated with QGSJet-II). 
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Figure 4. Muon production shower profiles from proton (red) 
and iron (blue) primaries, in {X',N') coordinates. The same 
showers used to build Fig. [T]are used here. 
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Figure 5. Average shower profiles for proton primaries at E = 
10" eV, with QGSJet-II, in (X',N') coordinates. Comparison 
between electromagnetic (in red) and muonic (in blue) shape fea- 
tures. The lines correspond to fits using a Gaisser-Hillas function 
(2 parameters). The fit results are given in the plot for the elec- 
tromagnetic (e.m.) and muonic (ji) profiles. 
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Figure 6. Correlation between the electromagnetic profile shape 
(characterized by R'"" ) and the muon production profile shape 
(represented by f). The results are shown for proton (red/full) 
and iron (blue/dashed) induced showers, at E = 10" eV, and with 
QGSJet-II. The shaded area shows the region within one sigma. 



half mean width and the rise of the profile happens faster 
than for proton induced showers. The impact of different 
hadronic interaction models in the shape was found to be 
very small |4]. 

This quasi-universal shape is compared to the electro- 
magnetic shape in Fig. |5] The MPD profile has a steeper 
growth and is more asymmetric, with respect to the shower 
maximum. Both profiles were fitted with a Gaisser-Hillas 
function, written in terms of a Gaussian width L and a 
parameter R, which is related with the asymmetry of the 
shower with respect to the shower maximum ||9i, such that. 



N' 
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X' 
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We conclude that the MPD profile can be well de- 
scribed with the same function as used for the electro- 



magnetic profile. In fact, it is worth noting that the full 
profile description is better achieved for the MPD profile. 
This is due to the fact that the energy deposit profile (the 
profile that is seen by the fluorescence detector and is in- 
timately related with the electromagnetic shower compo- 
nent) contains in the end-tail a non-negligible contribution 
from muons that is not accounted by the Gaisser-Hillas. 
Nevertheless, it is interesting to note that both distribu- 
tions, having different development mechanisms, share the 
same kind of structure. 

Through the study of the shape parameters (R, L) we 
were able to conclude that, similarly to what happened 
with the electromagnetic profile, one of the parameters 
retains almost all information on the type of primary al- 
lowing the other one to be fixed. However, while in the 
e.m. case it was U "- that, being related to the shower en- 
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ergy, could be fixed, in the MPD profile it is W that shows 
more capabilities to distinguish between proton and iron 
induced showers. Note that W is proportional to the in- 
tegral of the MPD profile which attains information about 
the total number of muons produced during the shower de- 
velopment. 

Having a single shape variable to characterise each 
profile, W and , the correlation between the shape 
of the MPD and the electromagnetic profile can be easily 
drawn. The relation between the shapes of the two pro- 
files for each individual event can be seen in fig. |6] Here, 
{U'"-,R^^) were fixed to its corresponding average values 
of the combined proton and iron distributions. The cor- 
relation is rather strong in the most populated region of 
{R'' '"-,W). Moreover, it is almost independent of the pri- 
mary mass composition, making it very useful for hybrid 
analysis. Indeed, whenever one of the profiles is measured 
accurately a prediction of the other profile shape can be 
established allowing, in this way to perform consistency 
tests to the shower description. 



4 Constraining IHadronic Interaction 
l\/lodels 

Finally, in this section we want to point out some strate- 
gies on how to disentangle the mass composition from 
the hadronic interaction models uncertainty. For this, it 
is useful to introduce the two-dimensional plots where the 
two first momenta of Xf^'^^ distribution, namely (^,^"^) and 
RMS (X^na';), are plotted for each mass composition sce- 
nario at a fixed energy [lol IJ. |. In these kind of plots each 
point represents a possible scenario and by considering the 
possible primary mass combinations we are, in a sense, 
probing the phase space for a given hadronic interaction 
model. 

Here, we want to extend this concept and combine 
two observables, that are sensitive to both the hadronic in- 
teraction models and primary mass composition: the av- 
erage total number of muons produced in the showeiQ, 
^A'^^, and the average electromagnetic shower maximum 
depth, (^f„"^)- The resulting plot is shown in Fig. Qfor 
three distinct hadronic interaction models: QGSJet-II.03 
(in black); EPOS 1.99 (in red) and SIBYLL2.1 Q (in 
green). Let us now evaluate the results of a single model, 
namely EPOS. The most extreme point at the left repre- 
sents a scenario with pure iron as mass composition, while 
the point at the right is for a scenario where there are only 
protons. The pure helium and nitrogen scenario are also 
marked for the case of EPOS. The remaining red dots are 
the results for bimodal mass composition scenario from 
one pure element to another Inside of this contour are 
all the remaining complex mass composition combination 
amongst the four species: p. He, N, Fe. 
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''this quantity is obtained performing the integral of the MPD profile 
up to the maximum. In this way, we obtain a quantity which is propor- 
tional to the total number of muons in the shower without any ground 
effects. For a fixed zenith angle the number of muons at ground gives 
similar results. 
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Figure 7. Average number of muons versus average electromag- 
netic X,';™; for E = 10" eV showers. Each single point in the plot 
represents a primary mass composition scenario. The results are 
shown for three distinct hadronic interaction models: QGSJet- 
II.03 (in black); EPOS1.99 (in red) and SIBYLL2.1 (in green). 
For EPOS the pure mass composition scenarios are marked with 
black crosses and the label below identifies the specie. 



The most interesting feature of this plot is that while 
the problem is degenerate when {Nf,'^ and are in- 

terpreted alone, here the hadronic interaction models can 
be distinguished independently of any mass composition 
scenario, provided, of course, that the experimental resolu- 
tions are better than the separation between models. More- 
over, these kind of analyses are important to constrain the 
allowed phase space for hadronic interaction models, in- 
creasing in this way our knowledge about the shower 

5 Conclusions 

The muon production longitudinal profile of air showers 
gives new primary mass composition variables which are 
fairly independent of the high energy hadronic interaction 
model: the X^„„^ and the shape variable W. The normali- 
sation of the profile gives also access to the total number 
of produced muons, which is known to be an important 
variable, both for primary composition and high energy 
hadronic interaction model studies. 

The MPD profile can be described using the param- 
eterisation that is used for the electromagnetic profile. 
Moreover, its shape displays an universal behaviour when 
expressed in terms of X' and A^'. 

Combining the information of the muonic longitudi- 
nal profile at production with the electromagnetic profile 
will give rise to extra variables, sensitive to the shower de- 
velopment characteristics. These will lead to a more pre- 
cise understanding of the shower development. In fact, 
we have shown that already some constrains on hadronic 
interaction models can be applied, independently of any 
mass composition scenarios, through the combination of 
the average electromagnetic shower maximum depth and 
the average number of muons produced in the shower. 
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With the higher number of available observables, using 
profiles which are independent of the detection conditions 
and more directly related to the hadronic cascade, cosmic 
rays become increasingly useful for the study of particle 
physics at the highest energies. 
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